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Abstract. We present several models for IRS 13E, an in- 
frared, mm and X-ray source in the Galactic Center. Our 
favored interpretation is that of an early-type binary with 
strong colliding winds emission. This naturally explains 
the observed X-ray count rate and the strong IR emis- 
sion lines, and has a distinct advantage over competing 
hypotheses based upon a single star or BH system. 



1. Introduction 

It is probable that Sgr A*, the compact, nonthermal ra- 
dio source at the Galactic Center (GC) is a 2-3x10 6 Mq 
black hole (see, e.g., Ghez et al. 1998). Deep in the po- 
tential well of Sgr A* and pervading the central parsec 
of the Milky Way, there exists a cluster of Hel and late- 
type stars (see, e.g. Sellgren et al. 1990 and Genzel et al. 
1996). In order to understand the evolution and dynamics 
of this unique stellar cluster, it is necessary to investigate 
individual stellar sources in detail. IRS 13, identified as 
a Pa-a, [Felll], Hel, and Hell line source (Stolovy et al. 
1999; Lutz, Krabbe, & Genzel 1993; Libonate et al. 1995; 
Krabbe et al. 1995), is a compact HII region, dominated 
by the source IRS 13E. 

Motivated by the Chandra results we speculate on the 
nature of IRS 13 and conclude that it is most likely an 
early-type binary system. This interpretation also sup- 
ports the other observations. 

2. Source Identification 

Due to the ~ 30 magnitudes of visual extinction towards 
the GC (Wade et al 1987; Blum, Sellgren, & Depoy, 1996), 
the only observations of IRS 13 are at radio, infrared and 
X-ray wavelengths. At A = 2 cm, the IRS 13 complex 
is one of the brightest sources in the central 10" of the 
Galaxy (Yusef-Zadeh, Roberts, & Biretta 1997). Obser- 
vations at 7mm and 13mm, where interstellar scatter- 
broadening is less significant, resolve IRS 13 into two 
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sources. Infrared K-band photometry also resolves IRS 
13E into two sources (IRS 13E1 and IRS 13E2), the former 
being identified with a Wolf-Rayet star of type WN9-10 
(Najarro et al. 1997, hereafter N97). The two sources are 
separated by 0.16" ~ 1000 AU in projection (assuming 
a GC distance of 8 kpc) with positional uncertainties of 
~ 0.1" (Ott, Eckart & Genzel 1999). This compares to 
a separation of ~ 500 AU for the mm sources although 
the positional uncertainties are consistent with these mm 
sources^] being the same as the infrared sources identi- 
fied by Ott, Eckart & Genzel (1999). Fig. |l] presents a 
schematic of the region, showing the mini-spiral and var- 
ious wind sources in the central 10 arcseconds. 

A recent Chandra observation of the GC (Baganoff et 
al. 1999) shows a compact spherical X-ray source appar- 
ently coincident with IRS 13, located ~ 3" ESE from Sgr 
A*. With ~ 100 counts after a 50 ksec observation, it is 
nearly as bright as Sgr A*. The only other point source 
in the central ~ 10" is a more asymmetrical object, with 
no obvious radio or infrared counterpart, that is located 
- 9" NNE of Sgr A*. 

3. Radio, sub-mm and IR observations 

(Zhao & Goss 1999) determined a spectral index of 0.9 ± 
0.2 for IRS 13E1, similar to that expected from a fully ion- 
ized stellar wind (Wright & Barlow 1975; Nugis, Crowther, 
& Willis 1998). Conversely IRS 13E2 has a flat mm spec- 
trum which may be nebular (Whiteoak 1994). Such a mix- 
ture of compact and extended emission is similar to the 
Luminous Blue Variable (LBV) r\ Car (see, e.g., Cox et al 
1995), which may have a binary companion (Damineli et 
al. 2000). 

A major feature of IRS 13 is the strength of some of 
its IR emission lines. IRS 13 has stronger 2.19/im Hell 
(N97) and Pa-a (Stolovy et al. 1999) lines than its nearby 
IRS 16 brethren and is one of the strongest Br-7 sources 
in the central parsec of the Galaxy (Libonate et al. 1995; 

1 Zhao & Goss (1999) label these source as IRS 13E and IRS 
13W; here, we will identify them as IRS 13E1 and IRS 13E2, 
respectively. 
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Fig. 1. A schematic of the central half parsec of the 
Galaxy, showing the IRS sources which are discussed in 
this paper as well as other early-type stars in the region. 
Shown in Crosshatch are some of the ionized gas compo- 
nents of the central mini-spiral. The location of Sgr A* is 
shown by an asterisk. The locations of IRS 13E1 and IRS 
13E2 are from Ott et al. (1999) and assume that IRS 7 is 
offset from Sgr A* by -0.1 and -5.69 arcseconds in x and 
y, respectively. 

Roberts, Yusef-Zadeh, & Goss 1996; Herbst et al. 1993). 
However, IRS 13 has a rather weak H92a line; Roberts 
et al (1996) thus suggest that the H92a source is hotter, 
denser and more compact than the Br7 source (Roberts, 
Yusef-Zadeh, & Goss 1996). Although its [Fell] emission 
is weak (Stolovy et al. 1999), IRS 13 is located at the local 
peak of [FcIII] emission (Lutz, Krabbe, & Genzel 1993), 
suggesting a local ionizing source. It should be noted, how- 
ever, that for some of the above observations, line blending 
(e.g. HeI(4-3) with Pa-a and HI(12-4) with the 1.644 ^m 
[Fell] line) can be a significant problem. 

From observations of the K-band Hel and Br7 emis- 
sion lines, N97 constructed non-LTE atmospheric models 
of IRS 13E1 and determined the parameters given in Ta- 
ble |l|. At first glance this mass loss rate is substantially 
more than the canonical upper limit of ~ 1 x 10~ 4 Mq 
yr _1 for WR stars (Nugis 1999). However, a number of 
arguments exist against such a high estimate. First N97 
did not include wind clumping or line-blanketing and as 
a result probably over estimated the mass loss rate by at 
least a factor of 2 (Morris et al. 2000). Second, assuming 
the parameters in Table |[ the classical stellar wind the- 



Table 1. Computed Model for IRS 13E1 (from N97) 



M (lO^MQyr- 1 ) 7.9 

(10 3 km s" 1 ) 1.0 

L, (1O 6 L ) 2.3 

T off (10 4 K) 2.9 

He/H > 500 



ory of Wright & Barlow (1975) predicts a 13 mm flux of 
~ 30 mJy, considerably higher than the observed ~ 10 
mJy for the combined IRS 13E complex, again suggest- 
ing an overestimate of M/woo by at least a factor of 2. 
Third, if the other GC wind sources have similarly over- 
estimated mass loss rates, hydrodynamical accretion sim- 
ulations (e.g. Coker & Melia 1997), which predict a large 
mass accretion rate onto Sgr A*, are brought more into 
agreement with other accretion models (Coker & Melia, 
2000; Narayan et al. 1998), which require a much lower 
accretion rate. For completeness, we also note that if the 
metallicity of IRS 13E is more than twice solar (see, e.g., 
Shields & Ferland 1994; Simpson et al. 1995), the maxi- 
mum possible WR mass loss rate could be higher (Maeder 
& Maynet 1987). 

If we assume the metallicity of the GC is twice so- 
lar then based on the estimated bolometric luminosity of 
IRS 13E1 the present mass and zero-age main-sequence 
(ZAMS) mass of IRS 13E1 are ~ 60 M Q and more than 
120 Mq, respectively (Schaller et al 1992). However, up- 
dated evolutionary models (e.g., Mowlavi et al 1998) show 
that a high metallicity star with ZAMS mass of £ 60 Mq 
will shed so much mass while on the main sequence that it 
would never reach the WR stage. The canonical assump- 
tion that the metallicity of the GC is twice solar, how- 
ever, is based on gas phase abundances. Carr, Sellgren, & 
Balachandran (2000) further argue that GC stellar abun- 
dances may in turn be overestimated due to higher stellar 
rotation rates. For example, high-resolution near-infrared 
spectra of IRS 7, a red supergiant located in the central 
parsec, show solar metallicity. 

4. What is IRS 13E? 

It is certain that IRS 13E is a luminous object embed- 
ded in a hot, dense medium, but more than this is not 
clear. However, it is possible to explore a few options which 
best fit the observations. First, the estimated luminosity 
and mass loss rate of IRS 13E1 are too large for a sin- 
gle evolved WR star; they are more like that of an LBV. 
In addition Libonate et al. (1995) noted that the K-band 
spectrum of IRS 13 resembles that of the LBV, P Cyg. 
Similarly, the FWHM of the Br-7 line (~ 350 km sec -1 ) 
as well as the extent of the emission region ( 0.12 pc) 
(Herbst et al. 1993) is consistent with the idea that IRS 
13E has recently undergone an LBV-like mass ejection, 
and also ties in well with the radio and sub-mm obser- 
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vations. However, in spite of this agreement its temper- 
ature and wind velocity are larger than a typical LBV. 
As has been suggested for IRS 16NE, IRS 16C, and IRS 
16SW (N97), it could therefore be a transition object just 
coming out of the LBV phase. This would be consistent 
with the v and M profiles of Garcia-Segura, Mac Low & 
Langer (1996). Humphreys & Davidson (1994) point out 
that a WN9/Ofpe star appears very much like an LBV 
at minimum brightness and thus the distinction between 
the two types is often blurred. However, the large He/H 
of IRS 13E1, compared to a He/H of - 1 for the IRS 16 
sources (N97), makes an LBV/WNL determination prob- 
lematic. Observational and theoretical counterarguments 
include the fact that the identity of WR 122, the calibrat- 
ing WNL source used by N97, has been called into doubt 
(Crowther & Smith 1999) and recent work (e.g., Langer et 
al. 1999) suggests mixing in massive stars is more efficient 
than previously thought, resulting in a larger He/H at the 
start of the WNL phase. Hence an LBV/WNL classifica- 
tion is certainly within reason (although see below). 

Second, we note that there are ^ 27 massive Hel 
stars in the GC (Blum, Ramirez, & Sellgrcn 1999). 
Given the observed frequency, /, of WR binary systems 
(12% £ / £ 50%; van der Hucht et al. 1981), it is likely 
that some of the GC Hel stars are binaries as well. In fact, 
IRS 16SW is thought to be an eclipsing binary with a pe- 
riod of - 10 days (Ott, Eckart & Genzel 1999). Thus it 
is possible that IRS 13E is also a binary system, contain- 
ing for example a WN10 primary with a ZAMS mass of 
~ 100 Mq and a somewhat less massive companion that 
is either an O star or another WR star, or, possibly, a 
massive compact object. 

5. Consequences 

Of the single and binary scenarios, the latter is preferred. 
As previously mentioned, the model of N97 underpredicts 
the K-band Hell emission for IRS 13E1 by a factor of 
~ 3; the colliding winds of a binary will produce more 
He + (Marchenko et al. 1997), potentially explaining this 
deficiency. Also, the ionizing flux from a massive O star 
would explain why IRS 13E stands out in Pa-a and [Felll] . 

An early-type binary system will also have strong 
shocks as a result of colliding stellar winds. The X-ray 
luminosity of a binary system will be brighter than that 
from a solitary star so that in order of increasing X-ray lu- 
minosity one qualitatively has (with everything else equal) 
WR, O, O/O, WR/O, WR/WR. However, this is modi- 
fied by the binary separation: if the binaries are too close, 
absorption will suppress the observed X-ray emission but 
if they are too far apart the shocks are largely adiabatic 
and do not produce significant additional X-ray luminos- 
ity (Pittard & Stevens 1997). 

The various WR sources of IRS 16 (including IRS 
16SW) do not appear as point sources in the Chandra 
observations. A solitary WR or even a WR with a B or 
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Fig. 2. Plot of a model WR/O spectrum detected by 
Chandra with a 50 ksec observation, assuming an in- 
trinsic (0.5-10.0 keV) X-ray luminosity of 0.5 Lq and a 
characteristic temperature of 1.5 keV (both of which are 
typical for colliding wind binaries, see ,e.g., Maeda, et 
al. 1999, Stevens, et al. 1996), and a column density of 
5.0 x 10 22 cm~ 2 . Chandra would see a total of ~ 100 pho- 
tons from this hypothetical binary system. 

later companion may not be visible with Chandra due 
to the large column density between here and the GC. 
Using XSPECf] and a Raymond-Smith thermal plasma 
model (which assumes optically thin X-ray line and con- 
tinuum emission; Raymond & Smith 1977) in ionization 
equilibrium one can simulate the Chandra spectrum of 
a solitary O-star placed at the GC. Assuming an ISM- 
corrected X-ray luminosity (0.5-10.0 keV) L x ~ 0.25 Lq 
(from L x /Lb i ~ 10~ 7 , see e.g. Waldron, et al. 1998), a 
characteristic temperature kT ~ 0.5 keV (representative 
of typical solitary O-stars, see Chlebowski, et al. 1989), 
an intervening column density Nh ~ 5 x 10 22 cm -2 (typ- 
ical for the GC, see e.g., Zylka et al. 1995), and solar 
abundances throughout, a 50 ksec Chandra observation, 
spanning the 0.5-10 keV band, would detect ^ 3 photons, 
and the O-star would not stand out above the background. 
Thus the IRS 16 sources blend into the diffuse background 
seen in the central ~ 5" of the image in Baganoff et al. 
(1999). Note that this hints that either the unseen IRS 
16SW companion does not have a significant stellar wind 
and thus is of type B or later (and therefore implies a 
lower system mass than estimated by Ott et al. 1999), or 
that the circumstellar absorption and/or binary separa- 
tion were unfavorable. 

In contrast, Fig. |^ shows a theoretical Chandra spec- 
trum for a hypothetical 0.5 Lq GC X-ray source similar 
to the WR/O binary j 2 Vel, which has a 78.5 day period. 
The number of photons detected by a 50 ksec observation 

2 Distributed and maintained by HEASARC 
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Fig. 3. Plot of a model WR/BH spectrum detected by 
Chandra with a 50 ksec observation, assuming an intrinsic 
luminosity of 0.1 Lq, a characteristic temperature of 5.0 
keV, and a column of 5.0 x 10 22 cm -2 . Chandra would see 
a total of about 45 photons from this hypothetical source. 

of such a source is ~ 100, consistent with the supposed 
IRS 13 source seen by Baganoff (1999). Note that the mm 
observations of IRS 13E2 probably do not correspond to 
the WR's companion but rather to the extended emission 
due to the the ejected nebula from the WR's LBV phase. 
Thus the binary separation is not constrained to match 
that of the resolved mm sources. 

Another possible interpretation of IRS 13 which would 
explain the Chandra observation follows from the sugges- 
tion by Gerhard (2000) that the GC Hel cluster is the 
remains of a disrupted cluster that formed tens of parsecs 
from Sgr A*. If true, a massive binary system may not 
have survived the infall. This raises the possibility that 
IRS 13E2 hides a compact object, such as a ~ 10 Mq 
black-hole, which is accreting the wind from IRS 13E1. 
That is, IRS 13E may be an extended X-ray system. This 
is particularly applicable in GC since dynamical friction 
over the lifetime of the Galaxy is likely to result in a high 
density of ~ 10 Mq black-holes in the central parsec (see, 
e.g., Morris 1993). One can estimate the X-ray luminosity 
of such a system using 



w 2 ( GM 

L x = eMc 2 - 

\ 2Dv 2 



(1) 



where e is the accretion efficiency (taken as 0.1), M is the 
mass of the black-hole, v is the relative velocity of the 
black-hole to the wind, and D is the distance between the 
WR star and the black-hole. If M — 10 Mq, v = 1000 km 
sec -1 , and D = 500 AU, then L x ~ 0.1 Lq. The spectrum, 
as shown in Fig. ||, would be shifted to slightly higher en- 
ergy, the flux would be less attenuated by the high column 
and thus the ~ 45 detected photons would be consistent 
with the IRS 13 detection of Baganoff (1999). Note that, 



if the ~ 10 mas yr _1 relative mm proper motions (Zhao 
& Goss 1999) are any indication, the system may be un- 
bound, making IRS 13E a unique transient X-ray source. 

In summary, whilst a single massive star of unusually 
high X-ray luminosity cannot be discounted, we believe a 
colliding winds binary system best fits the various obser- 
vations. A BH system also has some difficulties with recon- 
ciling the resolved mm sources (unless the second source is 
unassociated) . Long term monitoring at mm wavelengths 
as well as a long-integration Chandra observation would 
help determine precisely what type of object lurks in IRS 
13E. 
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